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Kerosene  is  widely  used  in  lamps  and  stoves  and  is  less  of  a  fire  hazard  than  other  fuels,  such  as  mineral 
spirits,  due  to  its  flash  point  being  above  40  C  and  low  vapor  pressure  at  normal  ambient  temperatures. 
At  room  temperature  and  pressure,  the  concentration  of  vapor  above  a  kerosene  spill  is  well  below  the 
lower  flammability  limit,  but  gasoline  addition  to  kerosene  will  increase  the  vapor  pressure  and  the 
ignitability  of  kerosene.  The  aim  of  this  paper  is  to  accurately  predict  the  evaporation  and  diffusion 
behavior  of  fuel  mixtures  of  gasoline  and  kerosene  in  an  arbitrary  mixture  ratio.  While  the  vapor 
pressures  of  gasoline  and  kerosene,  which  are  both  multi-component  liquid  fuels,  decrease  with  the 
progress  of  evaporation,  it  is  demonstrated  that  the  vapor  pressures  can  be  expressed  by  exponential 
functions  of  weight  loss  fractions.  The  vapor  pressures,  evaporation  rates  and  flash  points  of  gasoline/ 
kerosene  mixtures  can  be  modeled  accumrately  by  treating  them  as  two-component  mixtures,  using 
empirically-derived  constants  for  each  component,  which  this  paper  gives.  Furthermore,  a  predictive 
model  of  a  concentration  distribution  of  vapor  above  a  fuel  mixture  spill  on  a  floor  was  derived.  Ignition 
experiments  of  fuel  vapor  were  carried  out,  and  the  model  was  experimentally  verified.  The  method 
presented  in  this  paper  can  provide  useful  information  to  estimate  the  fire  hazard  in  scenarios  where  a 
fuel  mixture  of  gasoline  and  kerosene  in  an  arbitrary  mixture  ratio  is  spilt  on  a  floor. 

©  201 1  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

It  is  easy  to  obtain  kerosene  as  well  as  gasoline,  and  both  fuels 
are  often  used  as  an  accelerant  in  arson.  There  is  little  or  no 
difference  in  flame  heights  of  pool  fires  between  gasoline  and 
kerosene  [1],  On  the  other  hand,  kerosene  has  a  vapor  pressure 
much  lower  than  that  of  gasoline  and  cannot  form  a  flammable 
mixture  in  an  ambient  atmosphere  at  room  temperature.  In  some 
arson  cases,  mixtures  of  gasoline  and  kerosene  have  been  used  as 
accelerants.  Adding  kerosene  into  gasoline  decreases  the  vapor 
pressure  and  the  vaporization  rate,  reducing  the  size  of  the  area 
where  a  flammable  mixture  can  be  formed.  These  changes  are 
expected  to  be  very  sensitive  to  the  mixture  ratio  of  gasoline  and 
kerosene.  We  must  accurately  predict  where  a  flammable  mixture 
spreads,  in  order  to  estimate  the  fire  hazard  in  case  of  liquid  fuel 
spills  [2].  Shepherd  et  al.  discussed  the  change  in  flash  points  and 
vapor  pressures  of  kerosene  by  gasoline  contamination  [3],  Fingas 
empirically  demonstrated  that  most  crude  oil  and  petroleum 
products  evaporate  at  a  logarithmic  rate  with  respect  to  time  and 
presented  a  simple  model  for  predicting  the  weight  fraction  evapo¬ 
rated  as  a  function  of  temperature  [4-6],  However,  there  have  not 
been  any  previous  reports  of  quantitative  investigation  of  how  much 
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fire  hazard  decreases  by  adding  kerosene  into  gasoline  or  how  much 
fire  hazard  increases  by  adding  gasoline  into  kerosene. 

In  this  work,  fuel  mixtures  of  gasoline  and  kerosene  were  used 
as  samples,  and  the  change  in  evaporation  characteristics  of  these 
samples  by  evaporation  and  the  time  variation  of  a  concentration 
distribution  of  vapor  from  a  fuel  mixture  spill  on  a  floor  were 
measured.  “Weight  loss  fraction”  was  chosen  as  the  parameter 
signifying  progress  of  evaporation.  The  weight  loss  fraction  a  is 
given  by: 

a  =  (w0— w)/w0  (1) 

where  w0  is  the  initial  weight  of  fuel  mixture  (kg)  and  w  is  the 
current  weight  of  gasoline  (kg).  While  gasoline  and  kerosene  were 
multi-component  liquid  fuels,  they  were  assumed  to  be  a 
two-component  liquid  mixture  of  gasoline  and  kerosene,  from  which 
the  evaporation  model  of  fuel  mixtures  was  derived.  The  change  in 
concentration  of  fuel  mixtures  during  evaporation  was  predicted  by 
the  model,  and  compared  with  the  measured  values.  A  diffusion 
model  above  a  fuel  mixture  spill  on  a  floor  was  prepared  based  on 
the  unidimensional  diffusion  equation  [7],  and  a  method  for  predict¬ 
ing  a  concentration  distribution  of  flammable  vapor  generated  from 
a  fuel  mixture  of  gasoline  and  kerosene  in  an  arbitrary  mixture  ratio 
will  be  presented.  Furthermore,  ignition  experiments  of  fuel  vapor 
were  carried  out  and  the  test  results  were  compared  with  the 
ignition  prediction  in  order  to  verify  the  presented  method. 
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Nomenclature 

R 

gas  constant  (8.314  J/molK) 

S 

evaporative  area  (m2) 

a 

weight  loss  fraction  (-) 

t 

time  (s) 

ag 

weight  loss  fraction  of  gasoline  (-) 

U 

ignition  time  (s) 

ak 

weight  loss  fraction  of  kerosene  (-) 

T 

temperature  (K) 

C 

concentration  of  fuel  vapor  (vol%) 

Tf 

flash  point  (K) 

Qu 

concentration  of  fuel  vapor  on  the  surface  of  fuel 

V 

evaporation  rate  of  fuel  mixture  (kg/m2s) 

mixture  at  time  t  (vol%) 

vg 

evaporation  rate  of  gasoline  component  in  fuel 

Cz.t 

concentration  of  fuel  vapor  at  time  t  and  height 

mixture  (kg/m2  s) 

z  (vol%) 

Vk 

evaporation  rate  of  kerosene  component  in  fuel 

d 

thickness  of  fuel  mixture  spill  (mm) 

mixture  (kg/m2  s) 

D 

diffusion  constant  of  fuel  vapor  (m2/s) 

w 

current  weight  of  fuel 

hi 

ignition  height  (m) 

mixture  (kg) 

k 

mass  transfer  coefficient  (m/s) 

Wo 

initial  weight  of  fuel 

kg 

mass  transfer  coefficient  of  gasoline  (m/s) 

mixture  (kg) 

kk 

mass  transfer  coefficient  of  kerosene  (m/s) 

xg 

molar  fraction  of  gasoline  component  in  fuel 

m,  n 

evaporation  rate  constant  of  fuel  mixture  (-) 

mixture  (-) 

Mg 

molecular  weight  of  gasoline  (0.08  kg/mol) 

Xk 

molar  fraction  of  kerosene  component  in  fuel 

Mk 

molecular  weight  of  kerosene  (0.12  kg/mol) 

mixture  (-) 

Ng 

evaporative  molar  of  gasoline  in  fuel  vapor  (mol) 

Xgo 

initial  weight  fraction  of  gasoline  component  in  fuel 

Nk 

evaporative  molar  of  kerosene  in  fuel  vapor  (mol) 

mixture  (-) 

P 

total  vapor  pressure  of  fuel  mixture  (Pa) 

yg 

molar  fraction  of  gasoline  component  in  fuel  vapor  (-) 

P  0 

atmospheric  pressure  (101,300  Pa) 

yk 

molar  fraction  of  kerosene  component  in  fuel  vapor 

Pg 

partial  vapor  pressure  of  gasoline  component  in  fuel 

ocA,  0Cb,  pa,  pb  vapor  pressure  constant  (-) 

mixture  (Pa) 

aAg, 

aBg,  fAg,  fBg  vapor  pressure  constant  of  gasoline  (-) 

Pk 

partial  vapor  pressure  of  kerosene  component  in  fuel 

aAk,  «B,k.  pA,k<  Pbm  vapor  pressure  constant  of  kerosene  (-) 

mixture  (Pa) 

6 

standardized  time  (m  s/kg) 

2.  Experimental 

2.3.  Preparation  of  samples 

Motor  regular  gasoline  and  kerosene,  which  were  purchased  at 
a  gas  station  and  were  in  conformity  to  JIS  K  2202  and  2203,  were 
mixed  in  weight  ratios  of  two-to-one,  one-to-one  and  one-to-two, 
and  prepared  the  three  kinds  of  fuel  mixtures.  These  fuel  mixtures 
are  hereinafter  called  “2:1  fuel  mixture”,  “1:1  fuel  mixture”  and 
“1 :2  fuel  mixture",  respectively.  Six  degraded  fuel  mixtures  with 
different  weight  loss  fractions  from  0  to  0.25  in  increments  of  0.05 
were  prepared  by  leaving  the  fuel  mixtures  in  an  open  tray  (the 
base  area:  0.1  m2)  at  20  °C.  Vapor  pressures  and  flash  points  of 
seven  samples  including  a  virgin  fuel  mixture  of  each  mixture 
ratio  were  measured. 

2.2.  Vapor  pressure  measurement 

Vapor  pressures  of  the  degraded  fuel  mixtures  were  measured 
in  the  range  from  10  up  to  40  C  at  the  intervals  of  5  °C  using  an 
automated  vapor  pressure  tester  (Tanaka  Scientific  Ltd.  AVP-30D). 
A  sample  cylinder  was  filled  with  evaporated  sample  of  30  mL 
and  submerged  in  a  water  bath  maintained  at  each  measurement 
temperature.  The  cylinder  was  shaken  in  the  water  bath  for  8  min 
until  the  temperature  of  the  sample  inside  the  cylinder  was 
constant,  and  pressure  rise  in  the  cylinder  was  measured  as  a 
vapor  pressure. 

2.3.  Flash  point  measurement 

Flashpoints  were  measured  by  an  automated  closed  cup  flash 
point  tester  (Tanaka  Scientific  Ltd.  AVP-30D).  The  closed-cup  flash 
point  tests  were  carried  out  in  accordance  with  JIS  1<  2265-1. 
A  degraded  fuel  mixture  of  50  mL  was  put  into  a  sample  cup,  and 
the  cup  was  covered  with  a  lid.  A  pilot  flame  was  brought  close  to 


the  inner  space  of  the  cup  from  an  observation  hatch  at  the 
intervals  of  0.5  °C,  while  heating  the  sample  at  a  rising  rate  of 

1  °C/min.  The  flash  point  was  defined  as  the  lowest  temperature 
where  a  flash  was  observed.  The  procedures  after  covering  the 
cup  were  automated.  Each  sample  was  examined  three  times,  and 
the  lowest  temperature  was  considered  to  be  the  flash  point. 
Flash  points  of  less  than  —30  °C  cannot  be  measured,  because  a 
refrigerator  cannot  cool  liquid  refrigerant  to  less  than  -  40  C. 

2.4.  Evaporation  rate  measurement 

Evaporation  rates  were  measured  by  weight  loss  using  an 
electronic  balance.  The  balance  was  a  Sartorius  CP4202S,  capable 
of  measurements  to  0.01  g.  A  tarred  square  pan  (base  area:  0.1  m2) 
was  loaded  on  the  balance,  and  a  fuel  mixture  was  poured  into  the 
pan,  then  the  weight  loss  was  measured  with  time.  The  data  was 
recorded  by  PC  at  time  intervals  of  10  s  until  the  weight  loss 
fraction  reached  0.3.  The  experiments  were  conducted  in  a  fume 
hood.  The  fume  hood  fan  was  not  operated,  and  the  fuel  mixture 
was  evaporated  in  a  no  wind  condition.  The  amount  of  the  poured 
fuel  mixture  was  200  mL  (thickness  of  the  fuel  d:  2  mm).  A  K-type 
thermocouple  of  80  pm  diameter  was  put  on  the  base  of  the  pan, 
and  the  temperature  of  the  fuel  mixture  was  measured. 

2.5.  Evaporation  and  diffusion  experiment 

Evaporation  and  diffusion  experiments  were  carried  out  using 
a  cylindrical  vessel  of  60  cm  diameter  and  100  cm  height  with  an 
opening  on  the  top.  A  1:1  fuel  mixture  (thickness  of  fuel  spill  d:  1, 

2  and  3  mm)  was  gently  spilt  over  the  whole  surface  of  the  vessel 
base.  Concentrations  of  fuel  vapor  were  measured  on  various  time 
t  at  the  designated  height  h  (h  =  10,  20,  30,  40  and  50  cm).  Every 
experiment  was  carried  out  at  19  C.  The  experimental  apparatus 
of  evaporation  and  diffusion  from  a  fuel  spill  is  shown  in  Fig.  1. 
The  concentrations  of  fuel  vapor  were  measured  by  infrared 
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OUT 


Fig.  1.  Experimental  apparatus  of  evaporation  and  diffusion  from  gasoline  spill  on 
a  floor. 


hydro  carbon  detectors  (Riken  Keiki  Co.  Ltd.  Rl-415),  which  were 
calibrated  using  standard  iso-butane.  These  experiments  were 
carried  out  in  a  sufficiently  large  room,  and  we  can  consider  the 
concentration  of  fuel  vapor  outside  the  vessel  to  be  zero. 

2.6.  Ignition  experiment  of  fuel  vapor 

Ignition  experiments  were  carried  out  using  a  cubic  steel 
frame  vessel  of  1  m3  with  a  base  plate  in  order  to  verify  the 
concentration  distribution  of  fuel  vapor  predicted  by  the  evapora¬ 
tion  and  diffusion  model.  The  frame  was  covered  with  a  vinyl 
sheet.  A  1:1  fuel  mixture  of  1  L  (thickness  of  fuel  spill  d:  1  mm) 
was  gently  poured  into  the  vessel  using  vinyl  hose  and  spilt  over 
the  whole  surface  of  the  vessel  base.  Then,  a  spark  discharge  was 
generated  in  the  vessel  using  a  6000  V  neon  sign  transformer  (the 
gap  of  electrode  was  2  mm),  and  observed  the  existence  of 
ignition  to  the  fuel  vapor  evaporated  from  the  fuel  spill.  Sparks 
were  generated  at  various  ignition  heights  hi  in  a  fixed  time  fj 
(ti=5,  10,  20  and  30  mim)  after  the  spill.  The  heights  where 
sparks  were  generated  were  designed  at  10,  20,  30  and  40  cm.  The 
experimental  apparatus  is  shown  in  Fig.  2. 


3.  Results  and  discussion 

3.1.  Vapor  pressure  of  fuel  mixtures 

The  relation  between  vapor  pressure  p  [Pa]  and  temperature  T 
[1<]  in  a  multi-component  fuel  mixture  can  be  explained  by  the 
following  Clausius-Clapeyron  equation  [8,9]. 

In  p=A-B/T  (2) 

where  A,  B  are  the  Clausius-Clapeyron  constants  [-]  and  A  is 
equivalent  to  the  quotient  of  the  evaporation  enthalpy  by  the  gas 
constant.  While  the  Clausius-Clapeyron  equation  is  well  establish 
for  pure  substaces,  it  is  reported  that  the  equation  is  applied  for 
gasoline  and  the  constants  of  gasoline  are  linearly  related  to  the 


weight  loss  fraction  a  [10]. 

A  =  aAa+pA  (3) 

B  =  aBa+fB  (4) 

Substituting  Eqs.  (3)  and  (4)  into  the  Clausius-Clapeyron 
equation  yields 

In  p  =  (aA-ocB/T)a+([}A-/SB/T)  (5) 

p  =  exp(/?A-/?B/T)exp{(aA-aB/T)a}  (6) 


Clear  vinyl  sheet 

1  m 


Fig.  2.  Experimental  apparatus  to  ignite  the  fuel  vapor  above  the  mixed  fuel  spill 
on  a  floor. 


Table  1 

Evaporation  constants  used  in  the  prediction  model. 


Fuel 

ax  (-) 

M-) 

aB  (-) 

Hb(-) 

k  (m/s) 

Gasoline 

1.49 

21.0 

1.74  x  103 

3.03  x  103 

3.15  x  10~4 

Kerosene 

46.4 

28.7 

1.69  xlO4 

7.06  xlO3 

2.25  x  10-4 

0.6 


0.5 

0.4 

0.3 

0.2 

0.1 


> 


Weight  loss  fraction  a  [-] 


Fig.  3.  Variation  of  vapor  pressure  of  gasoline  and  kerosene  with  progress  of 
evaporation. 


Thus,  if  the  temperature  is  constant,  the  change  in  the  vapor 
pressure  by  evaporation  can  be  expressed  by  an  exponential 
function  of  the  weight  loss  fraction.  Since  kerosene  is  multi- 
component  hydrocarbon  mixture  as  well  as  gasoline,  it  was 
assumed  that  Clausius-Clapeyron  equation  was  established  in 
kerosene.  The  values  of  aA,  [SA,  aB  and  ftB  were  obtained  from  the 
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measured  vapor  pressure  of  degraded  samples  of  gasoline  and 
kerosene,  respectively.  These  values  are  called  “vapor  pressure 
constant”  and  shown  in  Table  1.  The  change  in  vapor  pressure  of 
gasoline  can  be  obtained  from  the  constants  using  Eq.  (6),  and 
change  curves  in  the  vapor  pressure  by  evaporation  are  shown  in 
Fig.  3.  After  the  following  section,  the  vapor  pressures  of  degraded 
fuel  mixtures  were  predicted  using  the  constants  presented  in 
Table  1. 

The  evaporation  model  of  a  multi-component  liquid  was 
presented  in  the  previous  work  [11].  While  gasoline  and  kerosene 
are  multi-component  fuel,  we  expressed  the  vapor  pressure  of 
each  fuel  by  an  exponential  function  of  a  weight  loss  fraction  and 
assumed  that  the  fuel  mixtures  are  two-component  liquids  of 
gasoline  and  kerosene.  Then,  the  variation  of  the  composition  of  a 
fuel  mixture  by  evaporation  was  predicted. 

The  composition  of  the  released  vapor  from  a  fuel  surface  was 
assumed  to  be  equal  to  that  of  the  gas  phase  that  was  in 
equilibrium  with  the  liquid  phase.  In  other  words,  the  molar 
fraction  of  each  component  in  vapor,  yg  and  yk  was  expressed  as 


the  proportion  of  each  partial  vapor  pressure  to  the  total  vapor 
pressure  as  shown  in  Eqs.  (7)  and  (8). 

Tg  =  xgPg/(xgPg  +  xkPkl  (?) 

yk  =  xkPk/(xgPg+xkPi<>  (8) 

where  xg  and  xk  are  the  molar  fraction  of  gasoline  and  kerosene 
component  in  liquid  fuel,  and  pg  and  pk  are  the  vapor  pressure  of 
gasoline  and  kerosene  component  in  the  pure  state. 

When  a  small  amount  of  a  fuel  mixture,  dw  [kg]  was  vapor¬ 
ized,  molar  amount  in  the  vapor  mixture  of  each  component  diVj 
[mol]  was  respectively  calculated  as  shown  in  Eqs.  (9)  and  (10). 


dNg  =ygdw/(ygMg+ykMk) 

(9) 

dNk  =  ykdw/(ygMg+ykMk) 

(10) 

where  Mg  is  the  molecular  weight  of  gasoline  (0.08  kg/mol)  and  Mk 
is  that  of  kerosene  (0.12  kg/mol).  Since  gasoline  and  kerosene  are 
multi-component  fuel,  the  average  molecular  weights  increase  with 


b 


Weight  loss  fraction  a  [-] 


— 

*k 

Fig.  4.  Variation  of  molar  fraction  of  each  component  in  liquid  phase  with  progress  of  evaporation  at  20  °C.  (a)  3:1  fuel  mixture,  (b)  2:1  fuel  mixture,  (c)  1 :1  fuel  mixture, 
(d)  1:2  fuel  mixture,  (e)  1:3  fuel  mixture. 
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the  progress  of  evaporation.  In  this  work,  the  molecular  weights 
were  assumed  to  be  constant  (0.08  and  0.12  kg/mol,  respectively)  in 
order  to  present  a  simple  evaporation  model. 

As  mentioned  in  Section  3.1,  the  vapor  pressure  of  gasoline 
and  kerosene,  pg  and  pk  were  given  by 

pg  =  exp(PAg-PBg/T)exp[(aAg-aBg/T)ag)  (1 1 ) 

Pk  =  exp(l3Ak-lfBk/T)exp[(ccAk-aBM/T)ak )  (12) 

where  ag  and  ak  are  the  weight  loss  fraction  of  gasoline  and 
kerosene,  respectively. 

The  variation  of  the  molar  fraction  of  gasoline  and  kerosene 
compartment  in  fuel  mixtures  with  the  progress  of  evaporation 
was  calculated  by  using  Eqs.  (7)— (12).  The  variations  of  the  molar 
fraction  in  3:1,  2:1,  1:1,  1:2  and  1:3  fuel  mixtures  at  20  °C  by 


evaporation  are  shown  in  Fig.  4.  The  variations  of  the  weight  loss 
fraction  of  gasoline  and  kerosene  are  shown  in  Fig.  5.  The  results 
showed  that  the  gasoline  component  were  preferentially  lost  in 
every  fuel  mixture. 

If  fuel  mixture  was  assumed  to  be  ideal  liquid,  then  vapor 
pressure  of  fuel  mixture  was  calculated  by  Eq.  (13). 

p=XgPg+XkP|< 

=  xg  exp(PAg-pBg/T)exp[(aAg-aBtg/T)ag) 

+xk  exp(fSAk-PBk/T)exp{(aAk-uB'k/T)ak }  (13) 

The  vapor  pressure  of  the  fuel  mixture  in  the  various  weight 
mixture  ratios  at  10-25  °C  was  calculated  by  using  Eq.  (13).  The 
calculated  and  measured  vapor  pressures  are  shown  in  Fig.  6.  In 
every  fuel  mixture,  the  calculated  vapor  pressure  agreed  well 
with  the  measured  one.  It  was  found  that  the  change  in  the  vapor 
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c  d 


— 

flg 

ak 

Fig.  5.  Variation  of  weight  loss  fraction  of  each  component  in  liquid  phase  with  progress  of  evaporation  at  20  °C.  (a)  3:1  fuel  mixture,  (b)  2:1  fuel  mixture,  (c)  1:1  fuel 
mixture,  (d)  1:2  fuel  mixture,  (e)  1:3  fuel  mixture. 
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pressure  with  the  progress  of  evaporation  of  fuel  mixture  could  be 
predicted  accurately  by  using  the  evaporation  model. 

Then,  the  changes  in  the  vapor  pressures  of  3:1,  2:1,  1:1,  1:2 
and  1 :3  fuel  mixtures  at  10-40  °C  by  evaporation  were  calculated. 
The  obtained  vapor  pressures  of  fuel  mixtures  were  plotted  by 
Clausius-Clapeyron  equation  as  shown  in  Fig.  7.  They  were  shown 
to  lie  almost  on  a  straight  line;  therefore,  the  vapor  pressure 
constants  of  the  fuel  mixtures  were  obtained  from  the  intercepts 
and  slopes  of  the  plots  and  were  shown  in  Table  2.  The  vapor 
pressure  of  the  degraded  fuel  mixture  at  an  arbitrary  temperature 
can  be  calculated  using  these  vapor  pressure  constants  and  Eq. 
(6). 

If  the  vapor  pressure  constant  of  gasoline  and  kerosene  are 
obtained,  vapor  pressure  of  fuel  mixture  in  an  arbitrary  mixture 
ratio  can  be  accurately  predicted  in  the  method  described  above. 


3.2.  Flash  point  of  fuel  mixtures 


At  the  flash  point,  the  concentration  of  vapor  on  a  fuel  surface 
is  equal  to  the  lower  explosive  limit  (LEL)  of  the  fuel  vapor-air 
mixture.  Since  vapor  generated  from  fuel  mixtures  is  almost 
composed  of  gasoline  component  as  shown  in  Fig.  5,  the  LEL  of 
a  fuel  vapor-  air  mixture  was  assumed  to  be  similar  to  that  of  a 
gasoline-  air  mixture, 1.4  vol%  [12].  The  predicted  flash  points  of 
degraded  samples  of  fuel  mixtures  were  obtained  by  calculating 
the  temperature,  where  the  concentration  of  vapor  on  a  fuel 
surface  is  equal  to  the  LEL,  i.e.,  the  vapor  pressure  is  equivalent  to 
1420  Pa,  using  Eq.  (6)  that  expresses  the  relation  between 
temperature  and  vapor  pressure.  The  flash  point  Tf  is  given  by 


j  _  _ Pb _ 

f—  aAa+pA- ln(1420) 


(14) 


The  change  in  the  flash  points  of  three  kinds  of  fuel  mixtures 
(2:1,  1:1  and  1:2  fuel  mixture)  by  evaporation  were  calculated  by 
using  Eq.  (14)  and  the  vapor  pressure  constants  presented  in 
Table  2.  The  measured  and  predicted  flash  points  are  shown 
in  Fig.  8.  The  predicted  flash  points  almost  agreed  with  the 
measured  results.  In  the  degraded  2:1  fuel  mixtures  of  a  <0.1, 
the  degraded  1:1  fuel  mixtures  of  a  <  0.05  and  the  1:2  fuel 
mixture,  the  flash  points  could  not  be  measured  because  they 
are  -30  °C  or  less.  The  flash  point  of  the  fuel  mixtures  rose  as 
evaporation  progressed,  and  the  flash  point  of  the  degraded  1 :2 
fuel  mixture  of  a  =  0.25  especially  rose  up  to  around  40  °C. 

The  vapor  pressure  of  fuel  mixture  in  an  arbitrary  mixture 
ratio  was  expressed  by  Eq.  (13).  Substituting  0  for  a  in  Eq.  (13) 
and  replacing  molar  fraction  of  gasoline  xg  with  initial  weight 
fraction  *go,  the  relation  between  the  flash  point  of  non-degraded 
fuel  mixture  Tf  and  initial  weight  fraction  Xg0,  is  given  by 


1.4  = 


Mi,X, 


kAgO 


+ 


MkXgo  +  M  g(  1  —  Xgo) 
JVJg(l-Xgo) 


MkXg0  +  ALg(l  — Xgo) 


exP(fiA,g-0B,g/Tf) 


exp  (/?/u-/WTf> 


05) 


The  change  of  the  flash  point  by  mixing  gasoline  into  kerosene 
was  calculated  by  using  Eq.  (15).  It  was  found  that  the  flash  point 
of  kerosene  rapidly  decreased  by  the  contamination  of  gasoline  as 
shown  in  Fig.  9. 


3.3.  Evaporation  rates  of  fuel  mixtures 

The  evaporation  rate  of  a  multi-component  liquid  is  in  propor¬ 
tion  to  the  vapor  pressure  p  [Pa]  as  follows  [13]: 

v=kpM/RT  (16) 


a  b 


Weight  loss  fraction  a  [-]  Weight  loss  fraction  a  [-] 


C 


0 - i - i -  - =C- - 

0  0.05  0.1  0.15  0.2  0.25  0.3 


■  Measured  ( T  =  25°C) 

a  Measured  ( T  =  20°C) 

♦  Measured  ( T  =  15°C) 

•  Measured  ( T  =  10°C) 

-  Calculated  (T  =  25°C) 

-  Calculated  (T  =  20°C) 

-  Calculated  (T  =  15°C) 

-  Calculated  (T  =  10°C) 


Weight  loss  fraction  a  [-] 


Fig.  6.  Measured  and  calculated  vapor  pressure  of  mixed  fuel,  (a)  2:1  fuel  mixture,  (b)  1:1  fuel  mixture,  (c)  1:2  fuel  mixture. 
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where  k  is  the  mass  transfer  coefficient  [m/s],  M  is  the  molecular 
weight  [kg/mol],  R  is  the  gas  constant  (8.314  J/mol  K)  and  T  is  the 
temperature  [K],  The  molecular  weight  of  gasoline  and  kerosene 
were  assumed  constant  at  0.08,  0.12  kg/mol,  respectively.  The 
mass  transfer  coefficients  of  gasoline  and  kerosene  were,  respec¬ 
tively,  obtained  from  the  slope  of  the  approximate  straight  lines 
by  plotting  the  evaporation  rate  against  the  value  of  pM/RT.  The 
mass  transfer  coefficients  were  shown  in  Table  1 . 

Correlations  of  weight  loss  fractions  with  evaporation  rates  of 
fuel  mixtures  were  shown  in  Fig.  10.  The  evaporation  rates  of  fuel 
mixtures  decreased  gradually  with  the  progress  of  evaporation 
because  more  volatile  components  were  lost  by  evaporation.  The 
evaporation  rate  of  a  multi-component  liquid  mixture  can  be 
expressed  by  the  sum  of  the  products  of  evaporation  rate  and 
molar  fraction  in  each  component  [14].  The  fuel  mixtures  of 
gasoline  and  kerosene  were  assumed  to  be  two-component 
liquids,  so  the  evaporation  rate  was  given  by 

v  =  xgVg  +xkuk  =  xgkgpgMg/RT +xkkkpkMk/RT  (17) 

where  kg,  kk  are  the  mass  transfer  coefficients  of  each  component 
[m/s].  The  coefficients  were,  respectively,  obtained  from  the 


measured  evaporation  rates  of  gasoline  and  kerosene  as  shown 
in  Table  1. 

The  changes  in  the  evaporation  rates  of  3:1,  2:1,  1:1,  1:2  and 
1:3  fuel  mixtures  by  evaporation  were  calculated  by  using 
Eq.  (17),  and  they  were  compared  with  the  measured  one  as 
shown  in  Fig.  10.  In  every  fuel  mixture,  the  calculated  evaporation 
rate  agreed  well  with  the  measured  one.  The  change  in  evapora¬ 
tion  rates  with  the  progress  of  evaporation  of  fuel  mixtures  could 
be  predicted  accurately,  assuming  the  fuel  mixtures  of  gasoline 
and  kerosene  were  the  two-component  liquids. 


3.4.  Amount  of  the  generated  vapor  from  a  fuel  mixture  spill 

The  changes  in  the  evaporation  rates  of  3:1,  2:1,  1:1,  1:2  and 
1:3  fuel  mixtures  at  10-25  °C  by  evaporation  were  calculated  by 
using  Eq.  (17),  and  presented  in  Fig.  11.  The  logarithm  of  the 
obtained  evaporation  rates  are  plotted  against  the  weight  loss 
fractions  in  Fig.  12.  The  plots  indicate  linear  relationships.  There¬ 
fore,  the  relation  between  a  weight  loss  fraction  and  an  evapora¬ 
tion  rate  can  be  approximately  expressed  by  an  exponential 


a  b 


1  IT  [K1]  1/7'  [K1] 


c  d 


l/T  [K1]  1/7’[K-'| 


l/T  [K1] 


a  =  0 

— A — 

a  =  0.05 

a  =  0. 1 

a  =  0.15 

a  =  0.2 

a  =  0.25 

Fig.  7.  Clausius-Clapeyron  plots  of  fuel  mixtures,  (a)  3:1  fuel  mixture,  (b)  2:1  fuel  mixture,  (c)  1:1  fuel  mixture,  (d)  1:2  fuel  mixture,  (e)  1:3  fuel  mixture. 
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function. 

v  =  mr"“  (18) 

where  m  and  n  are  the  evaporation  rate  constants.  These  con¬ 
stants  were  respectively  obtained  from  the  intercepts  and  slopes 
of  the  approximate  straight  lines  as  shown  in  Table  3. 

The  evaporation  rate  v  [kg/m2  s]  is  the  weight  loss  per  dt  and 
unit  area  and  is  expressed  as: 

v  =  -(l/S)(dw/dt)  (19) 

where  w  is  the  current  weight  of  a  fuel  mixture  [kg],  S  is  the 
evaporative  area  [m2]  and  t  is  time  [s],  Eq.  (19)  is  replaced  by 
Eq.  (18): 


-(l/S)(dw/dt)  =  mexp(-na) 

(20) 

where 

dw=  -w0da 

(21) 

Eq.  (20)  is  rewritten  to: 

exp(na)da  =  —  dt 
w0 

(22) 

The  boundary  condition  for  Eq.  (22)  is: 

a  =  0  at  t=  0 

(23) 

Table  2 

Vapor  pressure  constants  of  fuel  mixtures. 


Fuel 

ctA 

Pa 

a  b 

Pb 

3:1  fuel  mixture 

1.80 

20.8 

2.36  x  103 

3.03  x  103 

2:1  fuel  mixture 

2.00 

20.7 

2.69  x  103 

3.03  x  103 

1:1  fuel  mixture 

2.95 

20.5 

3.79  x  103 

3.03  x  103 

1:2  fuel  mixture 

8.76 

20.0 

7.04  x  103 

2.97  x  103 

1:3  fuel  mixture 

21.2 

19.2 

1.18  x  104 

2.86  x  103 

Eq.  (22)  can  be  solved  as: 

a  =  ^ln{mnSt/w0  +  l}  (24) 

Multiplying  the  both  sides  of  Eq.  (24)  by  w0,  the  amount  of 
evaporated  fuel  mixture  aw0  is  given  by: 

awo  =  ^ln{mnSt/w0-t-l}  (25) 

Thus,  the  amount  of  evaporated  fuel  mixture  is  expressed  by 
logarithm  of  time.  The  amount  of  evaporated  fuel  mixture  after  an 
arbitrary  time  passes  can  be  estimated  using  Eq.  (25). 


Initial  weight  fraction  of  gasoline  *goH 

Fig.  9.  Variation  of  flash  point  of  fuel  mixture  with  varying  initial  weight  fraction 
of  gasoline. 


a  b 
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■  Measured  flash  point 
-  Predicted  flash  point 


Fig.  8.  Measured  and  predicted  flash  points  of  degraded  gasoline,  (a)  2:1  fuel  mixture,  (b)  1:1  fuel  mixture,  (c)  1:2  fuel  mixture. 
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a  b 
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—  Measured  value 
-  Calculated  value 


Fig.  10.  Measured  and  calculated  evaporation  rate  from  a  fuel  spill  of  2  mm  thickness,  (a)  2:1  fuel  mixture,  (b)  1:1  fuel  mixture,  (c)  1:2  fuel  mixture. 
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Weight  loss  fraction  a  [-] 
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Weight  loss  fraction  a  [-] 


-  T  =  30°C 

-  T  =  25  °C 

-  T  =  20°C 

-  T  =  15°C 
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Fig.  11.  Evaporation  rate  curves  with  varying  temperature,  (a)  3:1  fuel  mixture,  (b)  2:1  fuel  mixture,  (c)  1:1  fuel  mixture,  (d)  1:2  fuel  mixture,  (e)  1:3  fuel  mixture. 
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Replacing  St/w0  by  a  standardized  time  9  [m2s/kg],  Eq.  (25)  is 
rewritten  by 

a  =  ^ln{mn0+l}  (26) 

By  illustrating  the  relation  between  a  standardized  time  6  and 
a,  an  evaporation  characteristics  curve,  which  is  independent  of 
the  evaporative  condition  such  as  the  spill  area  and  the  amount  of 
spilt  fuel  (i.e.,  the  thickness  of  spilt  fuel),  can  be  obtained.  The 
evaporation  characteristics  curves  of  3:1,  2:1, 1:1, 1:2  and  1:3  fuel 
mixtures  at  10-30  °C  are  shown  in  Fig.  13  using  m  and  n 
presented  in  Table  3.  The  amount  of  the  evaporated  fuel  mixture 
in  an  arbitrary  evaporative  condition  can  be  read  off  from  the 
curve  by  calculating  the  corresponding  9  value.  Thus,  if 
the  mixture  ratio  of  a  fuel  mixture  is  known,  the  amount  of  the 
generated  vapor  in  an  arbitrary  spill  condition  can  be  estimate  by 
the  method  described  above. 

3.5.  Evaporation  and  diffusion  experiment 

The  concentration  distribution  of  the  vapor  generated  from  a 
gasoline  spill  on  a  floor  can  be  predicted  by  the  unidimensional 
diffusion  model  [7].  The  vapor  pressures  of  fuel  mixtures  decrease 


with  the  progress  of  evaporation.  Therefore,  the  concentration  on 
the  surface  of  a  fuel  mixture  needs  to  be  changed  correspondingly 
with  the  progress  of  evaporation.  The  concentration  of  vapor  on 
the  surface  of  a  fuel  mixture  was  expressed  by  the  ratio  of  the 
vapor  pressure  to  the  atmospheric  pressure,  and  is  given  by 

P  exp(fA-fB/T)exp{(aA-ocB/T)a} 

Co'r  =  To  = - To -  (27) 

where  p0  is  the  atmospheric  pressure  (101,300  Pa),  a  is  the  weight 
loss  fraction,  which  is  obtained  from  the  accumulated  amount  of 
the  evaporated  fuel  mixture,  T  is  the  temperature  of  a  fuel 
mixture  [K]  and  xA,aB,fA  and  fB  are  the  vapor  pressure  constants 
of  fuel  mixtures,  the  constants  presented  in  Table  2. 

Let  us  consider  the  unidimensional  vertical  diffusion  in  the 
z-direction  through  an  infinite  plane  space,  of  thickness  A z  [m] 
(Fig.  14).  The  upward  volumetric  flux  at  time  t  and  height  z,  Jzt 
[m3/m2s  (  =  m/s)]  is  given  by 

Jz.t  =  D(Cz,t-Cz+Az,t)/Az  (28) 

where  D  is  the  diffusion  constant  [m2/s],  Czt  is  the  concentration 
of  vapor  at  time  t  and  height  z. 
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c  d 
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-  T  =  30°C 

-  r=25°C 

—  r  =  20°c 

—  r=i5°c 

r=  io°c 


Fig.  12.  a-\nv  plots  of  mixed  fuel,  (a)  3:1  fuel  mixture,  (b)  2:1  fuel  mixture,  (c)  1:1  fuel  mixture,  (d)  1:2  fuel  mixture,  (e)  1:3  fuel  mixture. 
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Taking  the  mass  balance  in  the  space  at  time  t  and  height  z,  the 
concentration  changes  of  vapor  A Czt  in  a  very  short  time  At  is 
expressed  as: 

A  Cz,t  =  (Jz_Azt-Jzl)At/Az  (29) 


where  Jz-Az,t  is  the  influx  at  height  z-Az,  Jz  t  is  the  outflux  at 
height  z. 

Substituting  Eq.  (28)  into  Eq.  (29)  gives 
A  CZit  =  D(Cz_hZt:—2Czlt-\-Cz+AZi[)At/ A z2  (30) 


Table  3 

Evaporation  rate  constants  of  fuel  mixtures. 


Temperature  T  (°C) 

Fuel 

3:1  fuel  mixture 

2:1  fuel  mixture 

1:1  fuel  mixture 

1 :2  fuel  mixture 

1 :3  fuel  mixture 

m  [-] 

n  [-] 

m  [-] 

n  [-1 

m  [-] 

n  I-] 

m  [-] 

n  [-] 

m  [-) 

n  [-] 

30 

5.02  x  10-4 

6.04 

4.63  x  10-4 

6.95 

3.73x10-“ 

9.59 

2.35x10-“ 

13.4 

1.43  x  10-“ 

15.1 

25 

4.32  x  10~4 

6.18 

3.99  x  10~4 

7.12 

3.23  x  10"4 

9.89 

2.08  x  10-“ 

14.1 

1.28x10-“ 

16.1 

20 

3.69  x  10~4 

6.32 

3.41  x  10~4 

7.29 

2.78  x  10"4 

10.2 

1.84x10-“ 

14.8 

1.15x10-“ 

17.3 

15 

3.14  x  10~4 

6.47 

2.91  x  10~4 

7.47 

2.38  x  10"4 

10.5 

1.62x10-“ 

15.6 

1.03  x  10-“ 

18.5 

10 

2.66  xlO-4 

6.62 

2.46  xlO-4 

7.64 

2.03  x  10-4 

10.8 

1.42x10-“ 

16.4 

9.17  xlO-5 

19.7 

a 


b 


Standardized  time  0  [m2s/kg] 


C 


Standardized  time  9  [m2s/kg] 


d 


Standardized  time  9  [m2s/kg] 


e 


Standardized  time  6  [m2s/kg] 


—  r  =  30°c 

—  r  =  25°C 

—  r  =  20°c 

—  r=i5°c 
T=  10°C 


Fig.  13.  Evaporation  characteristics  curves,  (a)  3:1  fuel  mixture,  (b)  2:1  fuel  mixture,  (c)  1:1  fuel  mixture,  (d)  1:2  fuel  mixture,  (e)  1:3  fuel  mixture. 


58 


K.  Okamoto  et  al.  /  Fire  Safety  Journal  49  (2012)  47-61 


The  concentration  distributions  of  fuel  vapor  were  calculated  by 
numerically  integrating  Eq.  (30)  with  an  initial  condition  Co=0  at 
f=0,  in  increasing  At  of  1  s  and  A z  of  0.5  cm.  The  vapor  concentra¬ 
tion  on  the  fuel  surface,  which  was  the  boundary  condition,  was 
calculated  by  Eq.  (27)  using  the  ambient  temperature.  Since  vapor 
generated  from  fuel  mixtures  was  almost  composed  of  gasoline 
component  as  shown  in  Fig.  5,  the  physical  properties  of  fuel 
mixtures  can  be  expected  to  be  similar  to  that  of  gasoline.  The 
diffusion  constant  of  fuel  vapor  D  was  assumed  to  be 
1.0  x  10~5  m2/s  as  well  as  that  of  gasoline  vapor  [7].  The  calculated 
concentrations  at  every  10  cm  height  from  the  base  up  to  50  cm 
height  in  the  cylindrical  vessel  for  an  hour  after  the  spill  of  a  1:1 
fuel  mixture  are  compared  with  the  measured  results  as  shown  in 
Fig.  15.  The  concentration  of  the  vapor  increased  at  each  height  over 
time.  The  increases  were  greater  in  the  case  of  the  fuel  depth  at 
2  mm  than  1  mm.  The  reason  is  that  the  amount  of  low-boiling 
point  component  is  more,  and  the  chemical  component  composi¬ 
tion  changes  more  slowly  at  2  mm  than  at  1  mm.  The  calculated 
vapor  concentrations  agreed  well  with  the  measured  results.  It  was 
found  that  the  concentration  of  the  vapor  generated  from  a  fuel 
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Fig.  14.  Vertical  diffusion  model  of  fuel  vapor. 


mixture  spill  on  a  floor  can  be  predicted  by  using  the  presented 
diffusion  model. 

The  concentrations  of  vapor  generated  from  fuel  mixture  spills 
in  various  mixture  ratios  were  predicted  using  the  vapor  pressure 
constants  presented  in  Table  2.  The  time  variations  of  the 
concentrations  in  every  height  are  shown  in  Fig.  1 6.  In  the  higher 
weight  fraction  of  gasoline,  the  vapor  concentration  was  esti¬ 
mated  to  become  higher.  Thus,  evaporation  and  diffusion  beha¬ 
vior  of  fuel  mixtures  of  gasoline  and  kerosene  in  an  arbitrary 
mixture  ratio  can  be  predicted  by  the  evaporation  model  of  the 
fuel  mixture  and  the  unidimensional  diffusion  model. 

3.6.  Ignition  experiment  of  fuel  vapor 

Ignition  experiments  were  carried  out  in  order  to  verify  the 
predicted  concentration  distribution  of  fuel  vapor  in  the  evapora¬ 
tion  and  diffusion  experiments.  "Ignition”  or  “non-ignition”  to  the 
vapor  above  the  fuel  mixture  spills  on  the  base  of  the  frame  vessel 
was  observed  when  a  spark  discharge  was  generated  in  various 
elapsed  times  at  various  heights.  The  experiments  were  carried 
out  at  23-24  °C.  The  combustion  behavior  in  the  experiment 
where  the  flammable  mixture  was  ignited  at  20  cm  height  in 
30  min  after  the  fuel  spill  were  shown  in  Fig.  17.  The  flame 
spherically  propagated  in  the  flammable  mixture  and  soon 
reached  the  upper  and  lower  end  of  the  flammable  mixture. 
Subsequently,  the  flame  was  no  longer  advanced  and  propagated 
only  horizontally  as  shown  in  Fig.  17.  From  the  results,  it  was 
confirmed  that  the  layer  concentration  distribution  with  a  nega¬ 
tive  gradient  was  formed  above  a  fuel  mixture  spill. 

The  vapor  concentrations  in  the  frame  vessel  were  calculated  by 
the  above-mentioned  evaporation  and  diffusion  model  of  fuel 
mixture.  Since  the  component  of  vapor  generated  from  a  fuel 
mixture  of  gasoline  and  kerosene  can  be  expected  to  be  similar  to 
that  of  gasoline  vapor,  the  flammable  range  of  the  fuel  vapor-air 
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Fig.  15.  Measured  and  predicted  concentration  of  vapor  when  1:1  fuel  mixture  was  spilt  on  a  floor,  (a)  d= 1  mm.  (b)  d= 2  mm.  (c)  d= 3  mm. 
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3:1  Fuel  mixture 
2: 1  Fuel  mixture 
1 : 1  Fuel  mixture 
1 :2  Fuel  mixture 
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Fig.  16.  Predicted  concentration  of  vapor  when  fuel  mixture  was  spilt  on  floor  with  1  mm  thickness  at  20  °C  (a)  z=  10  cm.  (b)  z=20  cm.  (c)  z=30  cm,  (d)  z=40  cm.  (e)  z=50  cm. 


Fig.  17.  Combustion  behavior  of  flammable  mixture  from  a  fuel  spill  on  a  floor  in  a  frame  vessel.  (hj=20  cm,  tj=30  min),  (a)  t=0  s.  (b)  t=  100  ms.  (c)  t=200  ms.  (d)  t= 300  ms. 
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mixture  was  assumed  to  be  from  1.4  to  7.6  vol%  as  well  as  that  of 
gasoline  vapor-air  mixture  [14],  and  the  possibility  of  ignition  was 
predicted.  The  evaporation  and  diffusion  model  were  verified  by  the 


Table  4 

Results  of  ignition  experiments  and  ignition  predictions. 


Ignition  height  hj  (cm) 

Ignition  time  U  (min) 

5 

10 

20 

30 

40 

X 

X 

X 

X 

X  (0%) 

X  (0%) 

X  (0.2%) 

X  (0.6%) 

30 

X 

X 

X 

X 

X  (0%) 

X  (0.1%) 

x  (0.9%) 

o  (1.8%) 

20 

X 

o 

o 

o 

X  (0.2%) 

o  (1.4%) 

o  (3.3%) 

o  (1.8%) 

10 

o 

X 

X 

X 

O  (4.0%) 

o  (6.5%) 

X  (7.9%) 

X  (8.1%) 

Upper  row:  results  of  ignition  experiments. 

Lower  row:  ignition  predictions  (predicted  concentration  of  fuel  vapor). 


UEL:  7.6vol% 

T.F.T.:  1.4vol% 

-  ti  =  5min 

-  tj  =  lOmin 

O:  Ignition 

x:  Non-ignition 

-  tj  =  20min 

-  tj  =  30min 

Fig.  18.  Results  of  Ignition  experiments  and  predicted  concentration  of  vapor  from 
mixed  fuel  spill  on  a  floor,  (a)  7—23  C.  (b)  7=24  C. 


comparison  of  the  results  of  ignition  experiments  with  the  ignition 
predictions.  The  experimental  results  and  ignition  predictions  are 
shown  in  Table  4,  and  the  predicted  vapor  concentrations  with  the 
experimental  results  are  shown  in  Fig.  18.  The  experimental  results 
at  20  and  40  cm  height  agreed  with  the  ignition  predictions  by  the 
calculated  concentration.  The  calculated  vapor  concentrations  at 
10  cm  height  in  10  min  and  at  30  cm  height  in  30  min  were  within 
the  flammable  range;  nevertheless,  ignition  was  not  observed  in  the 
corresponding  experiments.  The  calculated  vapor  concentrations  in 
these  conditions  were  close  to  the  flammable  limit,  therefore,  the 
experimental  errors,  such  as  the  liquid  temperature  change  by 
the  latent  heat  of  evaporation  and  the  turbulence  diffusion  when 
the  fuel  was  poured  on  the  vessel  base,  were  expected  to  cause  the 
difference  between  the  prediction  and  experimental  results.  While 
the  fuel  mixtures  were  multi-component,  the  molecular  weight, 
diffusion  constant  and  flammable  range  were  assumed  to  the 
specific  values.  It  was  supposed  that  the  assumption  was  also  a 
factor  that  caused  the  different. 

Thus,  it  was  confirmed  that  the  range  where  the  flammable 
mixture  is  formed  can  be  well  predicted  by  calculating  the  vapor 
concentration  distribution  of  fuel  mixtures  using  the  presented 
evaporation  and  diffusion  model.  The  fire  hazard  in  the  cases 
where  a  fuel  mixture  of  gasoline  and  kerosene  in  an  arbitrary 
mixture  ratio  is  spilt  on  a  floor  could  be  estimated  by  the 
presented  model. 


4.  Conclusion 

The  aim  of  this  paper  is  to  accurately  predict  the  evaporation 
and  diffusion  behavior  of  fuel  mixtures  of  gasoline  and  kerosene 
in  an  arbitrary  mixture  ratio.  The  evaporation  characteristics  of 
fuel  mixtures  and  concentration  distributions  of  vapor  from  a  fuel 
mixture  spill  on  a  floor  were  measured.  The  evaporation  model 
assuming  fuel  mixture  as  two-component  liquid  of  gasoline  and 
kerosene  was  presented.  The  vapor  pressures,  the  flash  points  and 
the  evaporation  rates  of  fuel  mixtures  were  predicted  by  using  the 
model.  The  calculated  values  were  compared  with  the  measured 
values,  and  the  change  in  the  evaporation  characteristics  of  fuel 
mixtures  by  evaporation  was  investigated  quantitatively.  The 
calculated  values  of  vapor  pressures,  flash  points  and  evaporation 
rates  corresponded  to  the  measured  values;  therefore,  it  was 
confirmed  that  the  change  in  the  evaporation  characteristics  by 
evaporation  can  be  predicted  by  the  presented  model. 

The  diffusion  model  above  a  fuel  mixture  spill  on  a  floor  based 
on  the  unidimensional  diffusion  equation  [7]  was  prepared, 
and  the  method  for  predicting  a  concentration  distribution  of 
flammable  vapor  generated  from  a  fuel  mixture  of  gasoline  and 
kerosene  in  an  arbitrary  mixture  ratio  was  presented.  The 
calculated  vapor  concentrations  agreed  well  with  the  measured 
results.  Experiments  to  ignite  the  vapor  above  a  fuel  mixture  spill 
on  the  base  plate  of  the  frame  vessel  were  carried  out,  and  the 
existence  of  the  ignition  by  a  spark  discharge  generated  in  various 
elapsed  times  at  various  heights  was  observed.  The  experimental 
results  agreed  well  with  the  ignition  predictions  by  the  calculated 
concentration.  From  these  results,  it  was  confirmed  that  the  range 
where  the  flammable  mixture  is  formed  can  be  well  predicted 
using  the  presented  evaporation  and  diffusion  model. 

Gasoline  and  kerosene  are  multi-component  fuel  mixture  and 
the  physical  properties  differ  with  each  product.  In  order  to 
predict  evaporation  and  diffusion  behaviors  of  fuel  mixtures 
accurately,  it  is  most  desirable  to  get  the  parameters  of  evapora¬ 
tion  characteristics,  aA,aB,pA,  j)B  and  k  in  each  case.  However,  the 
evaporation  and  diffusion  behavior  of  vapor  from  a  fuel  mixture 
spill  on  a  floor  can  be  predicted  approximately  using  the  para¬ 
meters  presented  in  Table  1  because  the  vapor  pressure  of 
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gasoline  is  defined  within  44-78  kPa  at  38  °C  by  the  standard  of 
fuels  (JIS  K  2202).  The  method  presented  in  this  paper  can  provide 
useful  information  to  estimate  the  fire  hazard  in  the  cases  where 
a  fuel  mixture  of  gasoline  and  kerosene  in  an  arbitrary  mixture 
ratio  is  spilt  on  a  floor. 
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